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■ The equation of state (EOS) of dense matter plays an important role in the super- 
\ nova phenomenon, the structure of neutron stars, and in the mergers of compact 

Q ' objects (neutron stars and black holes). During the collapse phase of a supernova, 

O ■ the EOS at subnuclear densities controls the collapse rate, the amount of delep- 

tonization and thus the size of the collapsing core and the bounce density. Prop- 
erties of nuclear matter that are especially crucial are the symmetry energy and 
the nuclear specific heat. The nuclear incompressibility, and the supernuclear EOS, 
play supporting roles. In a similar way, although the maximum masses of neutron 
^ \ stars are entirely dependent upon the supernuclear EOS, other important structural 

■ aspects are more sensitive to the equation of state at nuclear densities. The radii, 
moments of inertia, and the relative binding energies of neutron stars are, in partic- 

^ \ ular, sensitive to the behavior of the nuclear symmetry energy. The dependence of 

' the radius of a neutron star on its mass is shown to critically influence the outcome 
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of the compact merger of two neutron stars or a neutron star with a small mass 
black hole. This latter topic is especially relevant to this volume, since it stems from 
research prompted by the tutoring of David Schramm a quarter century ago. 
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1 Introduction 



The equation of state (EOS) of dense matter plays an important role in the 
supernova phenomenon and in the structure and evolution of neutron stars. 
Matter in the collapsing core of a massive star at the end of its life is com- 
pressed from white dwarf-like densities of about 10^ g cm^'^ to two or three 
times the nuclear saturation density, about 3-10^'' g cm~'^ or = 0.16 baryons 
fm~^. The central densities of neutron stars may range up to 5-10 Ug. At den- 
sities around Ug and below matter may be regarded as a mixture of neutrons, 
protons, electrons and positrons, neutrinos and antineutrinos, and photons. At 
higher densities, additional constituents, such as hyperons, kaons, pions and 
quarks may be present, and there is no general consensus regarding the prop- 
erties of such ultradense matter. Fortunately for astrophysics, however, the 
supernova phenomenon and many aspects of neutron star structure may not 
depend upon ultradense matter, and this article will focus on the properties 
of matter at lower densities. 

The main problem is to establish the state of the nucleons, which may be 
either bound in nuclei or be essentially free in continuum states. Neither tem- 
peratures nor densities are large enough to excite degrees of freedom, such 
as hyperons, mesons or quarks. Electrons are rather weakly interacting and 
may be treated as an ideal Fermi gas: at densities above 10'' g cm~^, they are 
relativisitic. Because of their even weaker interactions, photons and neutrinos 
(when they are confined in matter) may also be treated as ideal gases. 

At low enough densities and temperatures, and provided the matter does not 
have too large a neutron excess, the relevant nuclei are stable in the labora- 
tory, and experimental information may be used directly. The so-called Saha 
equation may be used to determine their relative abundances. Under more ex- 
treme conditions, there are a number of important physical effects which must 
be taken into account. At higher densities, or at moderate temperatures, the 
neutron chemical potential increases to the extent that the density of nucleons 
outside nuclei can become large. It is then important to treat matter outside 
nuclei in a consistent fashion with that inside. These nucleons will modify the 
nuclear surface, decreasing the surface tension. At finite temperatures, nuclear 
excited states become populated, and these states can be included by treat- 
ing nuclei as warm drops of nuclear matter. At low temperatures, nucleons in 
nuclei are degenerate and Fermi-liquid theory is probably adequate for their 
description. However, near the critical temperature above which the dense 
phase of matter inside nuclei can no longer coexist with the lighter phase of 
matter outside nuclei, the equilibrium of the two phases of matter is crucial. 

The fact that at subnuclear densities the spacing between nuclei may be of 
the same order of magnitude as the nuclear size itself will lead to substan- 
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tial reductions in the nuclear Coulomb energy. Although finite temperature 
"plasma" effects will modify this, the zero-temperature Wigner-Seitz approxi- 
mation employed by Baym, Bethe & Pethick (|1]) is usually adequate. Near the 
nuclear saturation density, nuclear deformations must be dealt with, includ- 
ing the possibilities of "pasta-like" phases and matter turning "inside-out" 
{i.e., the dense nuclear matter envelopes a lighter, more neutron-rich, liquid). 
Finally, the translational energy of the nuclei may be important under some 
conditions. This energy is important in that it may substantially reduce the 
average size of the nuclear clusters. 

An acceptable way of bridging the regions of low density and temperature, 
in which the nuclei can be described in terms of a simple mass formula, and 
high densities and/or high temperatures in which the matter is a uniform bulk 
fluid, is to use a compressible liquid droplet model for nuclei in which the drop 
maintains thermal, mechanical, and chemical equilibrium with its surround- 
ings. This allows us to address both the phase equilibrium of nuclear matter, 
which ultimately determines the densities and temperatures in which nuclei 
are permitted, and the effects of an external nucleon fluid on the properties 
of nuclei. Such a model was originally developed by Lattimer et al. (0) and 
modifled by Lattimer & Swesty (|^). This work was a direct result of David 
Schramm's legendary ability to mesh research activities of various groups, in 
this case to pursue the problem of neutron star decompression. After the fact, 
the importance of this topic for supernovae became apparent. 



2 Nucleon Matter Properties 



The compressible liquid droplet model rests upon the important fact that in 
a many-body system the nucleon-nucleon interaction exhibits saturation. Em- 
pirically, the energy per particle of bulk nuclear matter reaches a minimum, 
about -16 MeV, at a density = 0.16 fm~^. Thus, close to n^, its den- 
sity dependence is approximately parabolic. The nucleon-nucleon interaction 
is optimized for equal numbers of neutrons and protons (symmetric matter), 
so a parabolic dependence on the neutron excess or proton fraction, x, can 
be assumed. About a third to a half of the energy change made by going 
to asymmetric matter is due to the nucleon kinetic energies, and to a good 
approximation, this varies as (1 — 2x)^ all the way to pure neutron matter 
(a; = 0). The x dependence of the potential terms in most theoretical models 
can also be well approximated by a quadratic dependence. Finally, since at 
low temperatures the nucleons remain degenerate, their temperature depen- 
dence to leading order is also quadratic. Therefore, for analytical purposes, the 
nucleon free energy per baryon can be approximated as fbuikii^, x), in MeV, as 
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where av{n) = (2m*//i^)(7r/12n)^/^. The expansion parameters, whose values 
are uncertain to varying degrees, are the incompressibihty, Ks = 190 — 250 
MeV, the skewness parameter K'^ = 1780 — 2380 MeV, the symmetry energy 
coefficient = S^^ns) = 25 — 36 MeV, and the bulk level density parame- 
ter, ay{ns,x = 1/2) ~ (l/15)(m*(ns, X = l/2)/m) MeV~"^, where m* is the 
effective mass of the nucleon. Values for ■m*{ns,x = 1/2) /m are in the range 
0.7 — 0.9. The general definition of the incompressibihty is A' = 9dP/dn = 
dd^n^dfbuik/dn) /dn, where P is the pressure, and Kg = Kijig, 1/2). It is worth- 
while noting that the symmetry energy and nucleon effective mass (which 
directly affects the matter's specific heat) are density dependent, but these 
dependencies are difficult to determine from experiments. The parameters, 
and their density dependences, characterize the nuclear force model and are 
essential to our understanding of astrophysical phenomena. 

The experimental determination of these parameters has come from compar- 
ison of the total masses and energies of giant resonances of laboratory nu- 
clei with theoretical predictions. Some of these comparisons are easily illus- 
trated with the compressible liquid droplet model. In this model, the nucleus 
is treated as uniform drop of nuclear matter with temperature T, density 
and proton fraction Xj. The nucleus will, in general, be surrounded by and be 
in equilibrium with a vapor of matter with density no and proton fraction Xo- 
At low ambient densities n and vanishing temperature, the outside vapor van- 
ishes. Even at zero temperature, if n is large enough, greater than the so-called 
neutron drip density Ud ~ 1.6 ■ 10~^ fm~^, the neutron chemical potential of 
the nucleus is positive and "free" neutrons exist outside the nucleus. At finite 
temperature, the external vapor consists of both neutrons and protons. In ad- 
dition, because of their high binding energy, a— particles will also be present. 
The total free energy density is the sum of the various components: 



Here, Fh and represent the free energy densities of the heavy nuclei and the 
outside vapor, respectively. The energy densities of the electrons and photons, 
Ff, and F^, are independent of the baryons and play no role in the equilibrium. 
For simplicity, we neglect the role of a-particles in the following discussion 
(although it is straightforward to include their effect (Q)). 

In the compressible liquid drop model, it is assumed that the nuclear energy 



F = Fh + Fo + F^ + F, + F. 



(2) 
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can be written as an expansion in A^^^ and (1 — 2xi)^: 



Fh — Uriilfhulk + fsurf + fcoul + f trans] , 



(3) 



where the /'s represent free energies per baryon due to the bulk, surface, 
Coulomb, and translation, respectively. The bulk energy, for example, is given 
by Eq. (1). The surface energy can be parametrized as 



where R is the nuclear charge radius, h{T) is a calculable function of temper- 
ature, (Jo is the surface tension of symmetric matter, and ag = (n^ /^GttY^^Ss 
where 5*^ is the surface symmetry energy coefficient from the traditional mass 
formula. In this simplified discussion, the infiuence of the neutron skin (^, 
which distinguishes the "drop model" from the "droplet model", is omitted. 
The Coulomb energy, in the Wigner-Seitz approximation is 



where D{u) = 1 — 1.5u^^^ + 0.5m and u is the fraction of the volume occupied 
by nuclei. If the fractional mass of matter outside the nuclei is small, u ~ n/n,. 

It is clear that additional parameters, 5*^ and another involving the tempera- 
ture dependence of h, exist in conjunction with those defining the expansions of 
the bulk energy. The temperature dependence is related to the matter's critical 
temperature Tc at which the surface disappears. It is straightforward to demon- 
strate from the thermodynamic relations defining Tc, namely dPbuik/dn = 
and d'^Pbuik/drP = 0, that Tc oc \^K^. Therefore, the specific heat to be asso- 
ciated with the surface energy will in general be proportional to T^~^ oc K~^. 
About half the total specific heat originates in the surface, so Ks infiuences 
the temperature for a given matter entropy, important during stellar collapse. 

The equilibrium between nuclei and their surroundings is determined by min- 
imizing F with respect to its internal variables, at fixed n, Yg, and T. This is 
described in more detail in Refs. (H; ||), and leads to equilibrium conditions 
involving the pressure and the baryon chemical potentials, as well as a condi- 
tion determining the nuclear size R. The latter is analogous to the one found 
by Baym, Bethe & Pethick (|1|) who equated the nuclear surface energy with 
twice the Coulomb energy. The relations in Eqs. (4) and (5) lead to 



f,urf = 47ri?V(xi,T) = ATiR^h{T)[ao - (t,(1 - 2x,f] , 
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Experimental limits to Kg., most importantly from RPA analyses of the breath- 
ing mode of the giant monopole resonance give Kg = 230 MeV. It is also 
possible to obtain values from the so-called scaling model developed from the 
compressibile liquid drop model. The finite-nucleus incompressibility is 

K{A, Z) = {M/h')R'El , (7) 



where M is the mass of the nucleus and E^r is the breathing-mode energy. 
K{A, Z) is commonly expanded as 

K{A, Z)=K, + K^^rfA-^'^ + K^iP + K,^rfiI^A~^'^ + KcZ^A-^'^ , (8) 



and then fit by least squares to the data for E^r. Here the asymmetry I = 
1 — 2Z/A. For a given assumed value of Kg, and taking Kgurfi = 0, Pearson (||) 
showed that experimental data gave 



i^c - 15.4 - 0.065fr, ± 2 MeV , ~ 230 - 3.2/^, ± 50 MeV . (9) 

With minimal assumptions regarding the form of the nuclear force, Pearson (^) 
demonstrated that values of Kg ranging from 200 MeV to more than 350 MeV 
could be consistent with experimental data. 

But the liquid drop model predicts other relations between the parameters: 
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Here E{Z,A) is the total energy of the nucleus, and is equivalent to Eq. (3). 
The second of these equations simply expresses the equilibrium between the 
nucleus and the surrounding vacuum, which implies that the pressure of the 
bulk matter inside the nucleus is balanced by the pressure due to the curvature 
of the surface and the Coulomb energy. It can then be shown that 



Kc = -(3eV5r,)[8 + 27n3/- , 
Ks^rf = 47rrlao[9ny:/ao + 22 + 54nf , 
Ksurfi = ^nrlag[9ny:/as + 22 + 54n^/- ,(n,)/ir,] , 

Kj = 9[nlS':{ng) - 2nX{ns) - 9ntS:K)/- . (11) 

Primes denote derivatives with respect to the density. From these relations. 
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and again assuming Kgurfi = 0, Pearson demonstrated that an interesting cor- 
relation between Kg and K'^, where K'^ = — 27n^/^'^;^(n5), could be obtained: 

K'^ = -0.0860irf + (28.37 ± 2.65)K, . (12) 

Assuming ~ 190 - 250 MeV, this suggests that K'^ = 1780 - 2380 MeV, a 
potential constraint. Alternatively, eliminating K'^, one finds 

Ks = 137 A - 2<oMnla"Jao ± 23.2 MeV . (13) 

The second derivative of the surface tension can be deduced from Hartree-Fock 
or Thomas-Fermi semi-infinite surface calculations. For example, if a parabolic 
form of fhuik is used, one finds 

nyjoo = -6 (14) 

leading to Kg = 295.5 ± 23.2 MeV. In general, the density dependence of Sy 
will decrease the magnitudes of Kg and a'^ from the above values. 

It is hoped current experimental work will tighten these constraints. A short- 
coming of the scaling model is that, to date, the surface symmetry energy term 
was neglected. This is not required, however, and further work is necessary to 
resolve this matter. 

Because the surface energy represents the energy difference between uniformly 
and realistically distributed nuclear material in a nucleus, the parameter Sg 
can be related to the density dependence of Sy{n) and to Kg. If f^uik is as- 
sumed to behave quadratically with density around Ug, this relation can be 
particularly simply expressed (^: 

Sg_ _ _3_ai/2 } y/x 

Sy ^2 To i 1 - X 

Here, Sy = Sy{ng), ai/2 = {dr / d\n.n)n^/2 is a measure of the thickness of the 
nuclear surface and Tq = (47^^/3)"^/^ = R/A^^^. If Sy{n) is linear, then the 
integral is 2; if Sv{n) oc n^/^, then the integral is 0.927. Since ai/2 will be 
sensitive to the value of Kg, we expect the value of Sg/Sy to be also. 

Experimentally, there are two major sources of information regarding the sym- 
metry energy parameters: nuclear masses and giant resonance energies. How- 
ever, because of the small excursions in A^^^ afforded by laboratory nuclei, 
each source provides only a correlation between Sg and Sy. For example, the 
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total symmetry energy in the liquid droplet model (now explicitly including 
the presence of the neutron skin, see Ref. @) is 



E 



sym 
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(16) 



Evaluating a = din Sg/dlia Sy near the "best-fit" values Ssq and S^q, one finds 

a~2 + 5,o< A>i/V'5so-6, (17) 

where < A >^/^ for the fitted nuclei is about 5. Thus, as the value of Sy is 
changed in the mass formula, the value of Sg must vary rapidly to compensate. 

An additional correlation between these parameters can be obtained from the 
fitting of isovector giant resonances, and this has the potential of breaking 
the degeneracy of and Ss, because it has a different slope (H). Lipparini 
& Stringari (|^ used a hydrodynamical model of the nucleus to derive the 
isovector resonance energy: 
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where m* is an effective nucleon mass. This relation results in a slightly less- 
steep correlation between Ss and S'„, 



a = 2/m* + (3/5)5^0 < A >^/^ /^,o - 4 - 5 . 



(19) 



Unfortunately the value of m* is an undetermined parameter and this slope 
is not very different from that obtained from fitting masses. Therefore, uncer- 
tainties in the model make a large difference to the crossing point of these two 
correlations. A strong theoretical attack, perhaps using further RPA analysis, 
together with more experiments to supplement the relatively meager amount 
of existing data, would be very useful. 



3 The Equation of State and the Collapse of Massive Stars 



Massive stars at the end of their lives are believed to consist of a white dwarf- 
like iron core of 1.2-1.6 Mq having low entropy (s < 1), surrounded by layers 
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of less processed material from shell nuclear burning. The effective Chan- 
drasekhar mass, the maximum mass the degenerate electron gas can support, 
is dictated by the entropy and the average lepton content, Y^, believed to be 
around 0.41-0.43. As mass is added to the core by shell Si-burning, the core 
eventually becomes unstable and collapses. 

During the collapse, the lepton content decreases due to net electron capture 
on nuclei and free protons. But when the core density approaches 10^^ g cm~^, 
the neutrinos can no longer escape from the core on the dynamical collapse 
time (^. After neutrinos become trapped, is frozen at a value of about 0.38- 
0.40, and the entropy is also thereafter fixed. The core continues to collapse 
until the rapidly increasing pressure reverses the collapse at a bounce density 
of a few times nuclear density. 

The immediate outcome of the shock generated by the bounce is also depen- 
dent upon Yl. First, the shock energy is determined by the net binding energy 
of the post-bounce core, and is proportional to Yg^°/'^ (^. Second, the shock 
is largely dissipated by the energy required to dissociate massive nuclei in the 
still-infalling matter of the original iron core outside the post-bounce core. 
The larger the Yl of the core, the larger its mass and the smaller this shell. 
Therefore, the progress of the shock is very sensitive to the value of Yl. 

The final value of Yl is controlled by weak interaction rates, and is strongly 
dependent upon the fraction of free protons, Xp, which is proportional to 
exp(/ip/T), and the phase space available for proton capture on nuclei, which 
is proportional to fie — A, where fi = fin — f^p- Both are sensitive to the proton 
fraction in nuclei {xi) and are largely controlled by Yl- In addition, the specific 
heat controls the temperature which has a direct influence upon the free proton 
abundance and the net electron capture rate. In spite of the intricate feedback, 
nuclear parameters relating chemical potentials to composition, especially Sy 
and S's, are obviously important. 

As an example, consider (x 
Eqs. (3)- (5), one has 

fi = ASy{l - 2xi) - 



Recall that as oc Ss- Although the bulk and Coulomb terms alone (Eq. 20 
with as = 0) imply that fi for a given Xj rises with increasing Sy, the proper 
inclusion of the surface symmetry energy gives rise to the opposite behavior. 
This is illustrated in Fig. 1. 

Uncertainties in nuclear parameters can thus be expected to have an influence 
upon the collapse of massive stars, for example, in the collapse rate, the final 



f'n — f^p = —n-i dFn/dxi. With the model of 
' 72ne'^Dy ^^ao - ct,(1 - 2xi){l - 6xi) 



dXiUi I ((To - crJl - 2xi)2)i/3 



(20) 
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Fig. 1. Comparison of function of Xj for various assumed values of 

both including and excluding the effects of the surface symmetry energy. 

trapped lepton fraction, and the radius at which the bounce-generated shock 
initially stalls. Swesty, Lattimer & Myra (|10|) investigated the effects upon 
stellar collapse of altering parameters in a fashion constrained by nuclear sys- 
tematics. They found that as long as the parameters permitted a neutron 
star maximum mass above the PSR1913+16 mass limit (1.44 M©), the shock 
generated by core bounce consistently stalls near 100 km, independently of 
the assumed in the range 180-375 MeV and 5*^ in the range 27-35 MeV. 
Ref. ([T0|) also found that the final trapped lepton fraction is also apparently 
independent of variations in both Kg and S^. These results are in contrast 
to earlier simulations which had used EOSs that could not support cold, cat- 
alyzed 1.4 Mq stars, or in which Ss was not varied consistently with 5*^. The 
strong feedback between the EOS, weak interactions, neutrino transport, and 
hydrodynamics is an example of Mazurek's Law. 

In fact, the only significant consequence of varying 5*^, involved the pre-bounce 
neutrino luminosities. Increasing increases the electron capture rate (pro- 
portional to fie — 1^ and therefore increases the z/g luminosity during collapse, 
as shown in Fig. 2. Nevertheless, the collapse rate also increases, so that neu- 
trino trapping occurs sooner and the final trapped lepton fraction does not 
change. It is possible that large neutrino detectors such as Super-Kamiokande 
or SNO may be able to observe an enhanced early rise in neutrino luminosity 
from nearby galactic supernovae. 
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Fig. 2. The neutrino luminosities during infall as a function of the bulk symmetry 
energy parameter. 



4 The Structure of Neutron Steirs 



The theoretical study of the structure of neutron stars is crucial if new obser- 
vations of masses and radii are to lead to effective constraints on the EOS of 
dense matter. This study becomes ever more important as laboratory studies 
may be on the verge of yielding evidence about the composition and stiffness 
of matter beyond n^. To date, several accurate mass determinations of neutron 
stars are available, and they all lie in a narrow range (1.25 — 1.44 M©). There 
is some speculation that the absence of neutron stars with masses above 1.5 
M0 implies that M^ax for neutron stars has approximately this value. How- 
ever, since fewer than 10 neutron stars have been weighed, and all these are 
in binaries, this conjecture is premature. Theoretical studies of dense matter 
indicate that considerable uncertainties exist in the high-density behavior of 
the EOS largely because of the poorly constrained many-body interactions. 
These uncertainties are reflected in a significant uncertainty in the maximum 
mass of a beta-stable neutron star, which ranges from 1.5-2.5 Mq. 

There is some theoretical support for a lower mass limit for neutron stars in 
the range 1.1 — 1.2 Mq. This follows from the facts that the collapsing core 
of a massive star is always greater than 1 M© and the minimum mass of a 
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protoneutron star with a low-entropy inner core of ~ 0.6 Mq and a high- 
entropy envelope is at least 1.1 Mq. Observations from the Earth of ther- 
mal radiation from neutron star surfaces could yield values of the quantity 
Roo = R/ \/^ ~ 2GM/Rc^, which results from redshifting the stars luminosity 
and temperature. M—R trajectories for representative EOSs (discussed below) 
are shown in Figure 3. It appears difficult to simultaneously have M > IMq 
and Roo < 12 km. Those pulsars with at least some suspected thermal ra- 
diation generically yield effective values of R^o so small that it is believed 
that the radiation originates from polar hot spots rather than from the sur- 
face as a whole. Other attempts to deduce a radius include analyses (|1^ of 
X-ray bursts from sources 4U 1705-44 and 4U 1820-30 which implied rather 
small values, 9.5 < R^o < 14 km. However, the modeling of the photospheric 
expansion and touchdown on the neutron star surface requires a model de- 
pendent relationship between the color and effective temperatures, rendering 
these estimates uncertain. Absorption lines in X-ray spectra have also been 
investigated with a view to deducing the neutron star radius. Candidates for 
the matter producing the absorption lines are either the accreted matter from 
the companion star or the products of nuclear burning in the bursts. In the 
former case, the most plausible element is thought to be Fe, in which case 
the relation R ^ 3.2GM/c^, only slightly larger than the minimum possible 
value based upon causality, ( P^ ; |IBp is inferred. In the latter case, plausible 
candidates are Ti and Cr, and larger values of the radius would be obtained. 
In both cases, serious difficulties remain in interpreting the large line widths, 
of order 100-500 eV, in the 4.1 ± 0.1 keV line observed from many sources. 
A first attempt at using light curves and pulse fractions from pulsars to ex- 
plore the M — R relation suggested relatively large radii, of order 15 km (p!?]). 
However, this method, which assumed dipolar magnetic fields, was unable to 
satisfactorily reconcile the calculated magnitudes of the pulse fractions and 
the shapes of the light curves with observations. 

Prospects for a radius determination have improved in recent years, however, 
with the detection of a nearby neutron star, RX J185635-3754, in X-rays and 
optical radiation (|TB|). The observed X-rays, from the ROSAT satellite, are 
consistent with blackbody emission with an effective temperature of about 57 
eV and very little extinction. In addition, the fortuitous location of the star in 
the foreground of the R CrA molecular cloud limits the distance to D < 120 
pc. The fact that the source is not observable in radio and its lack of variabil- 
ity in X-rays implies that it is not a pulsar unlike other identified radio-silent 
isolated neutron stars. This gives the hope that the observed radiation is not 
contaminated with non-thermal emission as is the case for pulsars. The X-ray 
observations of RXJ185635-3754 alone yield Roo ~ 7.3(D/120 pc) km for a 
best-fit blackbody. Such a value is too small to be consistent with any neu- 
tron star with more than 1 Mq. But the optical fiux is about a factor of 2.5 
brighter than what is predicted for the X-ray blackbody, which is consistent 
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Fig. 3. M — R curves for the EOSs listed in Table 1. The diagonal lines represent 
two theoretical estimates (LP=Ref. ([TT|); RP=Ref. (|l2|)) of the locus of points for 
AI/I = 1.4% for extremal hmits of Pt, 0.25 and 0.65 MeV fm'^. The large dots 
on the M — R curves are the exact results. The region to the left of the contours 
labeled 0.65 is not allowed if current glitch models are correct (|l3|). 



with there being a heavy-element atmosphere (19). With such an atmosphere, 
it is found (|20|) that the effective temperature is reduced to approximately 50 
eV and -Roo is also increased, to a value of approximately 21.6(D/120 pc) km. 
Upcoming parallax measurements with the Hubble Space Telescope should 
permit a distance determination to about 10-15% accuracy. If X-ray spectral 
features are discovered with the planned Chandra and XMM space observa- 
tories, the composition of the neutron star atmosphere can be inferred, and 
the observed redshifts will yield independent mass and radius information. In 
this case, both the mass and radius of this star will be found. 

Furthermore, a proper motion of 0.34 " yr~^ has been detected, in a direction 
that is carrying the star away from the Upper Scorpius (USco) association ( pO]) . 
With an assumed distance of about 80 pc, the positions of RX J185635-3754 
and this association overlap about 800,000 years ago. The runaway OB star 
C Oph is also moving away from USco, appearing to have been ejected on 
the order of a million years ago. The superposition of these three objects is 
interesting, and one can speculate that this is not coincidental. If upcoming 
parallax measurements are consistent with a distance to RX J185635-3754 of 
about 80 pc, the evidence for this scenario will be strong, and a good age 
estimate will result. 
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Table 1 

Equations of state used in this work. Approach refers to the basic theoretical 
paradigm. Composition refers to strongly interacting components (n=neutron, 
p=proton, H=hyperon, K=kaon, Q=quark); all approaches include leptonic con- 
tributions. 



Symbol 


Reference 


Approach 


Composition 


FP 


m 


Variational 


np 


PS 


m 


Potential 


nvrO 


WFF(l-3) 




Variational 


np 


AP(l-4) 


(H) 


Variational 


np 


MS(l-3) 


m 


Field Theoretical 


np 


MPA(l-2) 




Dirac-Brueckner HF 


np 


ENG 




Dirac-Brueckner HF 


np 


PAL (1-6) 


m 


Schematic Potential 


np 


GM(l-3) 




Field Theoretical 


npH 


GS(l-2) 




Field Theoretical 


npK 


PCL(l-2) 


m 


Field Theoretical 


npHQ 


SQM(l-3) 


m 


Quark Matter 


Q iu,d,s) 



In this section, a striking empirical relationship is noted which connects the 
radii of neutron stars and the pressure of matter in the vicinity of Ug. In 
addition, a number of analytic, exact, solutions to the general relativistic TOY 
equation of hydrostatic equilibrium are explored that lead to several useful 
approximations for neutron star structure which directly correlate observables 
such as masses, radii, binding energies, and moments of inertia. The binding 
energy, of which more than 99% is carried off in neutrinos, will be revealed from 
future neutrino observations of supernovae. Moments of inertia are connected 
with glitches observed in the spin down of pulsars, and their observations yield 
some interesting conclusions about the distribution of the moment of inertia 
within the rotating neutron star. From such comparisons, it may become easier 
to draw conclusions about the dense matter EOS when firm observations of 
neutron star radii or moments of inertia become available to accompany the 
several known accurate mass determinations. 



4-1 Neutron Star Radii 

The composition of a neutron star chiefly depends on the nature of strong 
interactions, which are not well understood in dense matter. The several pos- 
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sible models investigated ([T5| ^ ) can be conveniently grouped into three broad 
categories: nonrelativistic potential models, field-theoretical models, and rel- 
ativistic Dirac-Brueckner-Hartree-Fock models. In each of these approaches, 
the presence of additional softening components such as hyperons, Bose con- 
densates or quark matter, can be incorporated. 



Figure 3 displays the mass-radius relation for several recent EOSs (the abbre- 
viations are explained in Table 1). Even a cursory glance indicates that in the 
mass range from 1 — 1.5 Mq it is usually the case that the radius has little 
dependence upon the mass. The lone exception is the model GSl, in which a 
kaon condensate, leading to considerable softening, appears. While it is gener- 
ally assumed that a stiff EOS leads to both a large maximum mass and a large 
radius, many counter examples exist. For example, MS3 has a relatively small 
maximum mass but has large radii compared to most other EOSs with larger 
maximum masses. Also, not all EOSs with extreme softening have small radii 
(viz., GS2). Nonetheless, for stars with mass greater than 1 Mq, only models 
with a large degree of softening can have Roc < 12 km. Should the radius of a 
neutron star ever be accurately determined to satisfy i?oo < 12 km, a strong 
case can be made for the existence of extreme softening. 
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12 13 
Radius (km) 

Fig. 4. Empirical relation between P and i? for various EOSs (see Table 1 for details). 
The upper and lower panels show results for gravitational masses of 1 and 1.4 
Mq, respectively. Symbols show PBr^^l^ in units of MeV fm~^ km~^/*^ at the three 
indicated fiducial densities. 
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It is relevant that a Newtonian polytrope with n = 1 has the property that 
the stellar radius is independent of both the mass and central density. In 
fact, numerical relativists have often approximated equations of state with 
n = 1 polytropes. An n = 1 polytrope has the property that the radius is 
proportional to the square root of the constant K in the polytropic pressure 
law P = Kp^~^^^"'. This suggests that there might be a quantitative relation 
between the radius and the pressure that does not depend upon the equation of 
state at the highest densities, which determines the overall softness or stiffness 
(and hence, the maximum mass). 

To make the relation between matter properties and the nominal neutron star 
radius definite. Fig. 4 shows the remarkable empirical correlation which exists 
between the radii of 1 and 1.4 M© stars and the matter's pressure evaluated 
at densities of 1, 1.5 and 2 n^. Table 1 explains the EOS symbols used in 
Fig. 4. Despite the relative insensitivity of radius to mass for a particular 
"normal" equation of state, the nominal radius Rm, which is defined as the 
radius at a particular mass M in solar units, still varies widely with the EOS 
employed. Up to ~ 5 km differences are seen in Ria, for example, in Fig. 4. 
This plot is restricted to EOSs which have maximum masses larger than about 
1.55 M0 and to those which do not have strong phase transitions (such as 
those due to a Bose condensate or quark matter). Such EOSs violate these 
correlations, especially for the case of 1.4 Mq stars. We emphasize that this 
correlation is valid only for cold, catalyzed neutron stars, i.e., it will not be 
valid for protoneutron stars which have finite entropies and might contain 
trapped neutrinos. The correlation has the form 



where P is the total pressure inclusive of leptonic contributions evaluated 
at the density n. An exponent of 1/4 was chosen for display in Fig. 4, but 
the correlation holds for a small range of exponents about this value. The 
correlation is marginally tighter for the baryon density n = 1.5ns and 2ns 
cases. Thus, instead of the power 1/2 that the Newtonian polytrope relations 
would predict, a power of approximately 1/4 is suggested when the effects of 
relativity are included. The value of the constant in Eq. (21) depends upon 
the chosen density, and can be obtained from Fig. 4. 

The exponent of 1/4 can be quantitatively understood by using a relativistic 
generahzation of the n = 1 polytrope, due to Buchdahl (pSf). For the EOS 



where is a constant, there is an analytic solution to Einstein's equations: 



R ~ constant • [-P(n)] 



0.23-0.26 



(21) 




(22) 
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e^'^gtt = {l-2l3){l-l3-u){l-(3 + v)-^ ; 

= g„ = (1 - 2/3)(l - /3 + - /3 - u)-\l - (3 + (3 cos Ar')-^ ; 
SirPG/c^ = A\\l -2(3){l-(3 + u)-^ ; 
87rpG/c^ = 2A^u{l -2(3){l-/3- 3'u/2)(l - (3 + u)'^ ; 

u = (3{Ar')-^ sin Ar' ; r = r'(l - /? + - 2/3)"^ ; 
A^ = 2887rp,Gc'-\l - 2/3)"^ R = 7i{l- /?)(! - 2P)-^A~\ (23) 



The free parameters of this solution are (3 = GM/Rc^ and the scale p*. Note 
that R (X p* ^^"^{1 + /9^/2 + ...), so for a given value of p*, the radius increases 
only very slowly with mass, exactly as expected from an n = 1 Newtonian 
polytrope. It is instructive to analyze the response of i? to a change of pressure 
at some fiducial density p, for a fixed mass M. One finds 



dlnR 
d\nP 



p,M 



d\nR 

dlnpt, 

1 + 



dlnP 



_£ 

dlnR 
din (3 



5 /P\ (l- /?)(!- 2/3) 

6 VpJ 2(1 -3/? + 3/32) 



(24) 



In the limit (3 ^ 0, P —>■ and ci In i^/ci In P 1/2, the value characteristic of 
an n = 1 Newtonian polytrope. Finite values of (3 and P render the exponent 
smaller than 1/2. If the stellar radius is about 15 km, p^, = 7r/(288i?^) ~ 
4.85 ■ 10"^ km-2. If the fiducial density is p ^ 1.5m;,n, ^ 2.02 ■ 10"^ km'^ 
(with nih the baryon mass), Eq. (22) implies that P ~ 8.5 ■ 10~^ km"^. For 
M = 1.4 Mq, the value of /? is 0.14, and d\nR/dlnP ~ 0.31. This resuh is 
mildly sensitive to the choices for p and R, and the Buchdahl solution is not a 
perfect representation of realistic EOSs; nevertheless, it provides a reasonable 
explanation of the correlation in Eq. (21). 

The existence of this correlation is significant because, in large part, the pres- 
sure of degenerate matter near the nuclear saturation density Ug is determined 
by the symmetry properties of the EOS. Thus, the measurement of a neutron 
star radius, if not so small as to indicate extreme softening, could provide an 
important clue to the symmetry properties of matter. In either case, valuable 
information is obtained. 

The specific energy of nuclear matter near the saturation density may be 
expressed as an expansion in the asymmetry (1 — 2x), as displayed in Eq. (1), 
that can be terminated after the quadratic term (P^). Leptonic contributions 
must be added to Eq. (1) to obtain the total energy and pressure; the electron 
energy per baryon is /g = (3 / 4:)hcx{3Tr'^nxY^^ . Matter in neutron stars is in 
beta equilibrium, i.e., /Xg — p^ + Pp = d{fbuik + fe)/dx = 0, so the electronic 
contributions may be eliminated to recast the pressure as (0) 
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1 ' \ 

K^ = 240 MeV 
S.,= 30 MeV 



1.5 



0.5 - 








Fig. 5. Left panel: M — R curves for selected PAL parametrizations (|2q ) showing the 
sensitivity to symmetry energy. The left panel shows variations arising from different 
choices of the symmetry energy at the nuclear saturation density = Sy{ns); the 
right panel shows variations arising from different choices of the density dependence 
of the potential part of the symmetry energy F{u) = Si;{n)/ S^ius) where u = n/ug. 



n 



S',{n){l - 2xy + 



xS„(n) 



n 



2x) + 



(- - 1) 



9ns ' hAug 
where x is now the beta equihbrium value. At the saturation density, 



(25) 



Ps = ns{l - 2xs)[nsS'^{ns){l - 2xs) + SyXs] 
where the equilibrium proton fraction at ng is 



(26) 



0.04 



(27) 



for Sv = 30 MeV. Due to the small value of Xg, one finds that Pg — n^S'^irig). If 
the pressure is evaluated at a larger density, other nuclear parameters besides 
Sv and S'yirig), become significant. For n = 2ng, one thus has 



P{2ng) ~ Ang[ngS',i2ng) + {Kg - K'j6)/9] 



(28) 
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If it is assumed that S'„(n) is linear in density, Kg ~ 220 MeV and K'^ ~ 2000 
MeV (as indicated in Eq. 12), the symmetry contribution is still about 70% 
of the total. 



The sensitivity of the radius to the symmetry energy is graphically shown by 
the parametrized EOS of PAL (|28D in Fig. 5. The symmetry energy function 
Sy{n) is a direct input in this parametrization. The figure shows the depen- 
dence of mass-radius trajectories as the quantities 5"^ and Sv{n) are alternately 
varied. Clearly, the density dependence of Sv{n) is more important in deter- 
mining the neutron star radius. Note also the weak sensitivity of the maximum 
neutron star mass to Sy. 



At present, experimental guidance concerning the density dependence of the 
symmetry energy is limited and mostly based upon the division of the nuclear 
symmetry energy between volume and surface contributions, as discussed in 
the previous section. Upcoming experiments involving heavy-ion collisions (at 
GSI, Darmstadt), which might sample densities up to ~ (3 — 4)ns, will be 
limited to analyzing properties of the symmetric nuclear matter EOS through 
a study of matter, momentum, and energy flow of nucleons. Thus, studies of 
heavy nuclei far off the neutron drip lines will be necessary in order to pin 
down the properties of the neutron-rich regimes encountered in neutron stars. 



4-2 Neutron Star Moments of Inertia and Binding Energies 



Besides the stellar radius, other global attributes of neutron stars are po- 
tentially observable, including the moment of inertia and the binding energy. 
These quantities depend primarily upon the ratio M/ R as opposed to details of 
the EOS, as can be readily seen by evaluating them using analytic solutions to 
Einstein's equations. Although over 100 analytic solutions to Einstein's equa- 
tions are known (p5|), nearly all of them are physically unrealistic. However, 
three analytic solutions are of particular interest in neutron star structure. 

The first is the well-known Schwarzschild interior solution for an incompress- 
ible fluid, p = pc, where p is the mass-energy density. This is mostly of interest 
because it determines the maximum compression jS = GM/ Rc^ for a neutron 
star, namely 4/9, based upon the pressure being finite. Two aspects of the 
incompressible fluid that are physically unrealistic, however, include the fact 
that the sound speed is everywhere infinite, and that the density does not 
vanish on the star's surface. 

The second analytic solution, Bl, due to Buchdahl (P^D, is described in Eq. (23). 
The third analytic solution (TolVII) was discovered by Tolman (^) in 1939, 
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and is the case when the mass-energy density p varies quadratically, that is, 

p = p,[l-(r/i?n. (29) 



In fact, this is an adequate representation, as displayed in Fig. 6 for neutron 
stars more massive than 1.2 Mq. The equations of state used are hsted in 
Table 1. The largest deviations from this general relation exist for models with 
extreme softening (GSl, GS2, PCL2) and which have relatively low maximum 
masses (see Fig. 3). It is significant that all models must, of course, approach 
this behavior at both extremes r ^ and r — >■ i?. 




.2 .4 .6 .8 .0 .2 .4 .6 .8 .0 .2 .4 .6 .8 .0 .2 .4 .6 .8 .0 .2 .4 .6 .8 .0 .2 .4 .6 .81.0 

r/R 



Fig. 6. Each panel shows mass-energy density profiles in the interiors of selected 
stars (masses indicated) ranging from about 1.2 Mq to the maximum mass (solid 
line) for the given equation of state (see Table 1). The thick black lines show the 
simple quadratic approximation 1 — {r/Kf'. 

Because the Tolman solution is often overlooked in the literature (for excep- 
tions, see, for example, Refs. (^; |H^)) it is summarized here. It is useful in 
establishing interesting and simple relations that are insensitive to the equa- 
tion of state. In terms of the variable x = r"^ / and the parameter (3, the 
assumption p = pc{l — x) results in pc = 15/?c^/(87rGi?^). The solution of 
Einstein's equations for this density distribution is: 
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1 - px{5 - 3x) 

^4 



^1 - 5f3/3)cos'^(f), 



iwi-w)/2 + C, (l)c 



f(5-3x) 



n 



pc^ + P cos (j) 
nihc'^ cos ( 



[X 



0) , C = tan-i v^/5/[3(l - 2/3)] , 
w = \og[x - 5/6 + Je'^/{3p)] , 1^;^ = w(x = 1) . (30) 



The central values of P/ pc^ and are 



P 



_2_ /3/Csc\2 /^y 



tan< 



tan< 



(31) 



This solution, like that of Buchdahl's, is scale-free, with the parameters (3 and 
Pc (or M or P). Here, n is the baryon density, rrib is the nucleon mass, and 
Csc is the sound speed at the star's center. When 0o = 7r/2, or /3 ^ 0.3862, 
Pc becomes infinite, and when f3 ~ 0.2698, Cgc becomes causal (i.e., c). Recall 
that for an incompressible fluid, Pc becomes infinite when (3 = 4/9. For the 
Buchdahl solution, Pc becomes infinite when (3 = 2/5 and the causal limit 
is reached when (3 = 1/6. For comparison, if causality is enforced at high 



densities, it has been empirically determined that (3 < 0.34 (|I5|; |T6|) 



The general applicability of these exact solutions can be gauged by analyz- 
ing the moment of inertia, which, for a star uniformly rotating with angular 
velocity fl, is 



R 



(87r/3) / r\p + P/c^)e^^-'^/\io/n)dr. 



(32) 



The metric function a; is a solution of the equation 

rf[r^e-(^+")/2^']/rfr + Ar^Lude-^^+^^/ydr = (33) 



with the surface boundary condition 

R , 



LJR = Q 



-UJ 



R 



2GI- 



(34) 



The second equality in the above follows from the definition of I and the TOV 
equation. Writing j = exp[— (i/ + A)/2], the TOV equation becomes 

j' = -AiTGr{P/c^ + p)3e^/c^ . (35) 
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Then, one has 




0.00 0.05 0.10 0.15 0.20 

M/R (Mo/km) 

Fig. 7. The moment of inertia / in units of MR'^ for the equations of state hsted in 
Table 1. line, Ibi, Ivii and Irp are approximations described in the text. 



Unfortunately, an analytic representation of oo or the moment of inertia for any 
of the three exact solutions is not available. However, approximations which 
are valid to within 0.5% are 

Iinc/MR^ ~ 2(1 - 0.87/3 - 0.3f3^)-^/5 , (37) 
Ibi/MR^ ~ (2/3 - 4/7r2)(l - 1.81/3 + 0.47/3^)-^ , (38) 
Itvii/MR^ ^2{1 - 1.1/3-0.6/3^)^77- (39) 

In each case, the small /3 limit reduces to the corresponding Newtonian results. 
Fig. 7 indicates that the Tolman approximation is rather good. Ravenhall & 
Pethick (|1^) suggested that the expression 

Irp/MR^ ~ 0.21/(1 - 2m) (40) 

was a good approximation for the moment of inertia; however, we find that 
this expression is not a good overall fit, as shown in Fig. 7. For low-mass 
stars (/3 < 0.12), none of these approximations is suitable, but it is unlikely 



22 



that any neutron stars are this rarefied. It should be noted that the Tolman 
approximation does not adequately approximate some EOSs, especially ones 
that are relatively soft, such as GM3, GSl, GS2, PAL6 and PCL2. 

The binding energy formally represents the energy gained by assembling 
baryons. If the baryon mass is rrib, the binding energy is simply BE = Nmb — 
M in mass units. However, the quantity rrih has various interpretations in 
the hterature. Some authors assume it is about 940 MeV/c^, the same as the 
neutron or proton mass. Others assume it is about 930 MeV/c^, corresponding 
to the mass of C^^/12 or Fe^^/56. The latter would yield the energy released in 
a supernova explosion, which consists of the energy released by the collapse of a 
white-dwarf-like iron core, which itself is considerably bound. The difference, 
10 MeV per baryon, corresponds to a shift of 10/940 ~ 0.01 in the value 
of BE/M. In any case, the binding energy is directly observable from the 
detection of neutrinos from a supernova event; indeed, it would be the most 
precisely determined aspect. 




0.5 1.0 1.5 2.0 2.5 3.0 

M (Me) 

Fig. 8. The binding energy of neutron stars as a function of stellar mass for the 
equations of state listed in Table 1. The predictions of Eq. (41) are shown by the 
shaded region. 

Lattimer & Yahil (^) suggested that the binding energy could be approxi- 
mated as 

BE ^ 1.5 ■ 10^^(M/Mq)2 ergs = O.O84(M/M0)2 M© . (41) 
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This formula, in general, is accurate to aboTit ±20%. The largest deviations 
are for extremely soft EOSs, as shown in Fig. 8. 




0.10 0.15 0.20 0.25 0.30 0.35 

GM/Rc^ 



Fig. 9. The binding energy per unit gravitational mass as a function of compactness 
for the equations of state listed in Table 1. The shaded region shows the prediction 
of Eq. (42) with ±5% errors. 

However, a more accurate representation of the binding energy is given by 

BE/M ~ 0.6/5/(1 - 0.5/3) , (42) 



which incorporates some radius dependence. Thus, the observation of super- 
nova neutrinos, and the estimate of the total radiated neutrino energy, will 
yield more accurate information about M/R than about M alone. 

In the cases of the incompressible fluid and the Buchdahl solution, analytic 
results for the binding energy can be found: 



B£WM=^(55^-;rr^)-i, (43) 

BEbi/M=(1 - 1.5/3)^1 - 2P(1 - 13)'^ - 1 . (44) 

The analytic results, the Tolman VII solution, and the fit of Eq. (42) are 
compared to some recent equations of state in Fig. 9. It can be seen that. 
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except for very soft cases like PS, PCL2, PAL6, GSl and GS2, both the Tolman 
VII and Buchdahl solutions are rather realistic. 



4-3 Crustal Fraction of the Moment of Inertia 

In the investigation of pulsar glitches, many models associate the glitch size 
with the fraction of the moment of inertia which resides in the star's crust, 
usually defined to be the region in which dripped neutrons coexist with nu- 
clei. The high-density crust boundary is set by the phase boundary between 
nuclei and uniform matter, where the pressure is Pt and the density n^. The 
low-density boundary is the neutron drip density, or for all practical purposes, 
simply the star's surface since the amount of mass between the neutron drip 
point and the surface is negligible. We define AR to be the distance between 
the points where the density is fit and zero. One can apply Eq. (32) to de- 
termine the moment of inertia of the crust alone with the assumptions that 
P/c^ « p, m{r) ~ M, and ujj ~ ujr in the crust. One finds 

AI ~ — / or e dr ~ — — / r dP , (45) 

3 n J ^ 3GM n J ' ^ ^ 

R-AR 

where M is the star's total mass and the TOV equation was used in the last 
step. In the crust, the fact that the EOS is of the approximate polytropic form 
P ~ Kp^/^ can be used to find an approximation for the integral / r^dP, viz. 

Pt 

j r^dP ~ PtR^ 



Since the approximation Eq. (42) gives / in terms of M and R. and cur/Q, is 
given in terms of I and R in Eq. (34), the quantity AI/I can thus be expressed 
as a function of M and R with the only dependence upon the equation of state 
(EOS) arising from the values of Pf and n^; there is no explicit dependence 
upon the higher-density EOS. However, the major dependence is upon the 
value of Pt, since rit enters only as a correction. We then find 

A7 _ 287TPtR^ (1 - 1.67/3 - 0.6(3^) 
~T ~ 3Mc^ P 

In general, the EOS parameter Pt, in the units of MeV fm^"^, varies over 
the range 0.25 < Pt < 0.65 for realistic EOSs. The determination of this 
parameter requires a calculation of the structure of matter containing nuclei 



1 + 



2Pt (l + 7/3)(l-2/3) 



(46) 



1 + 



2Pt (l + 7/3)(l-2/3) 



ritmiyC'^ 



/52 



.(47) 
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just below nuclear matter density that is consistent with the assumed nuclear 
matter EOS. Unfortunately, few such calculations have been performed. Like 
the fiducial pressure at and above nuclear density which appears in the relation 
Eq. (21), Pt should depend sensitively upon the behavior of the symmetry 
energy near nuclear density. 

Choosing nt = 0.07 fm~^, we compare Eq. (47) in Fig. 3 with full structural cal- 
culations. The agreement is good. We also note that Ravenhall & Pethick ( 1T2|) 
developed a different, but nearly equivalent, formula for the quantity AI/I as 
a function of M, R, Pf and /it, where fit is the neutron chemical potential at 
the core-crust phase boundary. This prediction is also displayed in Fig. 3. 

Link, Epstein & Lattimer ( |13D established a lower limit to the radii of neutron 
stars by using a constraint derived from pulsar glitches. They showed that 
glitches represent a self-regulating instability for which the star prepares over 
a waiting time. The angular momentum requirements of glitches in the Vela 
pulsar indicate that more than 0.014 of the star's moment of inertia drives 
these events. If glitches originate in the liquid of the inner crust, this means 
that AI/I > 0.014. A minimum radius can be found by combining this con- 
straint with the largest realistic value of Pt from any equation of state. Stellar 
models that are compatible with this constraint must fall to the right of the 
Pt = 0.65 MeV fm~^ contour in Fig. 3. This imposes a constraint upon the 
radius, namely that R > 3.6 -|- 3.9M/M0 km. 



5 The Merger of a Neutron Star with a Low-Mass Black Hole 

The general problem of the origin and evolution of systems containing a neu- 
tron star and a black hole was first detailed by Lattimer & Schramm (p9D, 
although the original motivation was due to Schramm. Although speculative 
at the time, Schramm insisted that this would prove to be an interesting topic 
from the points of view of nucleosynthesis and gamma-ray emission. The con- 
temporaneous discovery ( ^D] ) of the first-known binary system containing twin 
compact objects, PSR 1913+16, which was also found to have an orbit which 
would decay because of gravitational radiation within 10^° yr, bolstered his 
argument. Eventually, this topic formed the core of Lattimer's thesis (^T]), and 
the recent spate of activity, a quarter century later, in the investigation of the 
evolution and mergers of such compact systems has wonderfully demonstrated 
Schramm's prescience. 

Compact binaries form naturally as the result of evolution of massive stellar 
binaries. The estimated lower mass limit for supernovae (and neutron star 
or black hole production) is approximately 8 Mq. Observationally, the num- 
ber of binaries formed within a given logarithmic separation is approximately 
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constant, but the relative mass distributions are uncertain. There is some in- 
dication that the distribution in binary mass ratios might be flat. The number 
of possible progenitor systems can then be estimated. Most progenitor systems 
do not survive the more massive star becoming a supernova. In the absence 
of a kick velocity it is easily found that the loss of more than half of the mass 
from the system will unbind it. However, the fact that pulsars are observed to 
have mean velocites in excess of a few hundred km/ s implies that neutron stars 
are usually produced with large "kick" velocities originating in the supernova 
explosion. In the case that the kick velocity, which is thought to be randomly 
directed, opposes the star's orbital velocity, the chances of the post-supernova 
binary remaining intact increases. In addition, the separation in a surviving 
binary will be reduced significantly. Subsequent evolution then progresses to 
the supernova explosion of the companion. More of these systems survive be- 
cause in many cases the more massive component explodes. But the surviving 
systems should both have greatly reduced separations and orbits with high 
eccentricity. 

OrbI;al Decay Time Reduction 




e 



Fig. 10. The reduction of the gravitational radiation orbital decay time as a function 
of initial orbital eccentricity. The dashed line is the inverse of the Peters ( ^2|) / 
function; the dotted line shows f~^^^, which reasonably reproduces the exact result. 

Gravitational radiation then causes the binary's orbit to decay, such that cir- 
cular orbits of two masses Mi and M2 with initial semimajor axes a satisfying 

a < 2.8[MiM2(Mi + M2)/M|]i/^R0 , (48) 
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will fully decay within the age of the Universe (~ 10^° yr). Highly eccentric 
orbits will decay much faster, as shown in Fig. 10. The dashed curve shows 
the inverse of the factor (^2j) by which the gravitational wave luminosity of an 
eccentric system exceeds that of a circular system: 

/ = (1 + 73eV24 + 37eV96)(l - e^)-^/^^ ^^g^ 

Because the eccentricity also decays, the exact reduction factor is not as strong 
as 1//. A reasonable approximation to the exact result is f~'^^^, shown by the 
dotted line in Fig. 10. The coefficient 2.8 in Eq. (48) is increased by a factor 
of f~^^^^ or about 2 for moderate eccentricities. 

Ref. argued that mergers of neutron stars and black holes, and the subse- 
quent ejection of a few percent of the neutron star's mass, could easily account 
for all the r-process nuclei in the cosmos. Ref. ( ^91) is also the earliest reference 
to the idea that compact object binary mergers are associated with gamma-ray 
bursts. A later seminal contribution by Eichler, Livio, Piran & Schramm (143|) 
argued that mergers of neutron stars occur frequently enough to explain the 
origin of gamma-ray bursters. 

Since the timescale of gamma-ray bursts, being of order seconds to several min- 
utes, is much longer than the coalescence timescale of a binary merger (which 
is of order the orbital frequency at the last stable orbit, a few milliseconds), 
it is believed that a coalescence involves the formation of an accretion disc. 
Although neutrino emission from accreting material, resulting in neutrino- 
antineutrino annihilation along the rotational axis, has been proposed as a 
source of gamma rays, it seems more likely that amplification of magnetic 
fields within the disc might trigger observed bursts. In either case, the lifetime 
of the accretion disc is still problematic, if it is formed by the breakup of the 
neutron star near the Roche limit. Its lifetime would probably be only about 
a hundred times greater than the orbital frequency, or less than a second. 
However, this timescale would be considerably enhanced if the accretion disc 
could be formed at larger radii than the Roche limit. A possible mechanism is 
stable mass transfer from the neutron star to the black hole that would cause 
the neutron star to spiral away as it loses mass (|^; |45|) . 

The classical Roche limit is based upon an incompressible fluid of density p 
and mass M2 in orbit about a mass Mi. In Newtonian gravity, this limit is 

RRoche,Ne^t = (M, /OmOliT p)'/^ = 19.2{M, /M^p,,)'/^ km, (50) 

where pis = p/10^^ g cm~^. Using general relativity. Fishbone (^) found 
that at the last stable circular orbit (including the case when the black hole 
is rotating) the number 0.0901 in Eq. (50) becomes 0.0664. In geometrized 
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units, Rjioche/Mi = 13(14.4) (Mfpi5/M|)~^/^, where the numerical coefficient 
refers to the Newtonian (last stable orbit in GR) case. In other words, if the 
neutron star's mean density is pis = 1, the Roche limit is encountered beyond 
the last stable orbit if the black hole mass is less than about 5.9 Mq. Thus, for 
small enough black holes, mass overflow and transfer from the neutron star to 
the black hole could begin outside the last stable circular orbit. And, as now 
discussed, the mass transfer may proceed stably for some considerable time. In 
fact, the neutron star might move to 2-3 times the orbital radius where mass 
transfer began. This would provide a natural way to lengthen the lifetime of 
an accretion disc, by simply increasing its size. 

The final evolution of a compact binary is now discussed. Define q = irins/ Mbh-, 
/i = rrinsMBH/M, and M = Mbh + ^ns-, where m^^ and Mbh are the neutron 
star and black hole masses, respectively. The orbital angular momentum is 

J2 = Gii'^Ma = GM^aq^il + qf . (51) 

We can employ Paczyhski's ( |17| ) formula for the Roche radius of the secondary: 

i?,/a = 0.46[g/(l + g)]^/^ (52) 



or a better fit by Eggleton (|48|): 



Ri/a = 0.49[.6 + g-^/^ ln(l + q^^^)]-^ . (53) 



The orbital separation a at the moment of mass transfer is obtained by setting 
Re = R, the neutron star radius. For stable mass transfer, the star's radius 
has to increase more quickly than the Roche radius as mass is transferred. 
Thus, we must have, using Paczyriski's formula, 

dlnR dlnRp din a 1 , , 

= a> — = — + - (54) 



dlnTJins dlnnins dlnnins 3 

for stable mass transfer, a is defined in this expression, and is shown in Fig. 11 
for a typical EOS. If the mass transfer is conservative, than J = Jew, where 

. 32 G^^ ^I'^M'/^ _ 32 G'/^ q^M'/^ 

5 c5 a7/2 - 5 c5 (l + g)W/2 ^^^^ 



and 



j _ _^ qjl-q) 
J~2a^q{l + q)- 
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Fig. 11. din R/dlnnins = a (solid curve) and neutron star radius R (dashed curve) 
as functions of neutron star mass m„s for a typical dense matter equation of state. 



This leads to 



. a 5 \ 32 q^M^ 
^ '.2 6 ~ V - ~Y^(l + g)a4 



(57) 



Since mns < Mbh, ^7 < 0, and the condition for stable mass transfer is simply 
q < 5/6 + a/2. For moderate mass neutron stars, a ~ 0, so in this case 
the condition is simply q < 5/6, which might even be achievable in a binary 
neutron star system. Had we used the more exact formula of Eggleton, Eq. 
(53), we would have found q < 0.78. Note that it has often been assumed that 
R oc in such discussions (^5)) , which is equivalent to a = —1/3. This is 

unjustified, and results in the upper limit q = 2/3 which might inappropriately 
rule out stable mass transfer in the case of two neutron stars. 

A number of other conditions must hold for stable mass transfer to occur. 
First, the orbital separation a at the onset must exceed the last stable orbit 
around the black hole, so that a > 6GMbh/c^, or 



Second, the tidal bulge raised on the neutron star must stay outside of the 
black hole's Schwarzschild radius. Kochanek gives an estimate of the 
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height of the tidal bulge needed to achieve the required mass loss rate: 

nl/3 



Ar 



(59) 



where Pt is a dimensionless parameter of order 1 and = G^I'^M^I'^ j c?!'^ is the 
orbital frequency. For q we use the equality in Eq. (57), which is equivalent to 



R 



sh 



2GMbh/c^ < a 



R-Ar. 



(60) 



Finally, so that the assumption of a Roche geometry is valid, it should be pos- 
sible for tidal synchronization of the neutron star to be maintained. Bildstein 
& Cutler ( p9D considered this, and derived an upper limit for the separation 
at which tidal synchronization could occur by integrating the maximum 



"syn 



torque on the neutron star as it spirals in from infinity and finding where the 
neutron star spin frequency could first equal the orbital frequency. They find 



400M3 



(-)°. 



(61) 



which translates to 



400(- 



GMbh\5 a (l + qY 



Rc^ > R. 



< 1. 



(62) 



Next we consider the effect of putting some of the angular momentum into an 
accretion disc. Following the discussion of Ref. (|49|) , we assume an accretion 
disc contains an amount of angular momentum that grows at the rate 



j, = -(l-/)M='/V/^(l + g)-^g. 



(63) 



where / is a parameter, taken to be a fit to the numerical results of Hut & 



Paczyhski (|50D : 



/ = 5gi/V3 - 3g^/72 



2/3, 



(64) 



We then find the new condition for angular momentum conservation to be 



J + J, = J, 



GW ■ 



(65) 



which yields 



1 
6 



> 



32 q^M^ 
'T c^(l + g)a4 



(66) 
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Therefore, the new condition for stable mass transfer is 

(g'-/)/(i + g)<a/2-i/6. 



(67) 



The case / = 1 corresponds to neglecting the existance of an accretion disc. 

It remains to determine when an accretion disc is likely to form. Initially, 
matter flowing from the neutron star to the black hole through the inner 
Lagrangian point passes close to the black hole and falls in. However, as 
the neutron star spirals away, the accretion stream trajectory moves out- 
side the Schwarzschild radius. When the trajectory doesn't even penetrate the 
marginally stable orbit, an accretion disc will begin to form. Particle trajectory 



computations of the Roche geometry by Shore, Livio & van den Huevel ([5T|) 
suggest that its closest approach to the black hole is 

i^c = a(l + g)(0.5 - 0.227 In g)^ (68) 

Equating Rc to GGMbh/c^ yields 

(0.5 - 0.2271n#(l + g)>6^ff^^. (69) 



These constraints and allowed regions for stable mass transfer are shown in 
Fig. 12. Apparently, stable mass transfer ceases when rUns ~ 0.14 if the 
formation of an accretion disc is ignored. If the effects of disc formation are 
included, the stable mass transfer ceases when rrins ~ 0.22 Mq. In both cases, 
the neutron star mass remains above its minimum mass (about 0.09 Mq for 
the equation of state used here). Thus, the neutron star does not "explode" 
by reaching its minimum mass. 

Fig. 13 shows the time development of the orbital separation a and the neutron 
star's mass and radius during the inspiral and stable mass transfer phases. 
Solid lines are calculated assuming there is no accretion disc formed, while 
dashed lines show the effects of accretion disc formation. The time evolutions 
during stable mass transfer are obtained from Eq. (57) and Eq. (66), using 
rrins = qM/ (1 + g)^. With disc formation, the mass transfer is accelerated and 
the duration of the stable mass transfer phase is shortened considerably. Also, 
the neutron star spirals out to a smaller radius, and does not lose as much 
mass, as in the case when the accretion disc is ignored. 

Therefore, if stable mass transfer can take place, the timescale over which 
mass transfer occurs will be much longer than an orbital period, and lasts 
perhaps a few tenths of a second. This is not long enough to explain gamma- 
ray bursts. However, we have also seen the likelihood that an accretion disc 
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Fig. 12. The dark and light shaded regions show the binary masses for which mass 
transfer in a black hole-neutron star binary will be stable in the absence of, and the 
presence of, an accretion disc. The constraints Eq. (58) (a > 6Mbh)i Eq. (60) (Tidal 
bulge OK), Eq. (62) (Tidal synchronization), and Eq. (69) (Accretion disc forms) 
are shown by the appropriately labelled curves. The parallel, diagonal, dashed lines 
show evolutionary tracks for the labelled total BH+NS masses, beginning in each 
case with nins = 1.5 Mq. 

forms is quite large. Furthermore, the accretion disc extends to about 100 km. 
Even though this is considerably less than Ref. (^5]) estimated, the lifetime of 
such an extended disc is considerable. To order of magnitude, it is given by 
the viscous dissipation time, or 



aCgH 



(70) 



Here D is the radial size of the disc, a is the disc's viscosity parameter, Cg 
is the sound speed and H is the disc's thickness. Note that ~ QH where 
Q = 2'k/P = ^JgMbh / D'^ is the Kepler frequency. Thus, 

Since the magnitude of a is still undetermined, and usually quoted (0) to be 
about 0.01, and H is likely to be of order i?, we find Tmsc ~ 230 s for our case. 
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Fig. 13. The separation of a 1.5 Mq neutron star with a 3 Mq black hole during 
a merger is indicated by the dot-dashed line during inspiral and by a solid line in 
the outspiral during stable mass transfer. Other curves show the neutron star mass 
and radius during the stable mass transfer (outspiral) phase. Solid (dashed) lines 
are computed by ignoring (including) the effects of an accretion disc. 

This alleviates the timescale problem for these models. Numerical simulations 

of such events are in progress, and it remains to be seen if a viable gamma-ray 
burst model from neutron star-black hole coalescence is possible. If it is, a 
great deal of the credit should rest with Dave. 

We thank Ralph Wijers for discussions concerning accretion disks. 
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